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Abstract 

Quasi-elastic light scattering spectroscopy is regularly used to examine the dynamics of dilute 
solutions of diffusing mesoscopic probe particles in fluids. For probes in a simple liquid, the 
light scattering spectrum is a simple exponential; the field correlation function (q, r) of the 
scattering particles is related to their mean-square displacements X 2 = ((Ax(t)) 2 ) during r via 
g^(q,r) = exp(— ^q 2 X 2 ). However, demonstrations of this expression refer only to identical 
Brownian particles in simple liquids, and show that if the form is correct then it is also true for 
all r that g^(q,r) = exp(— Tt), a pure exponential in r. In general, g^'(q,r) is not a single 
exponential in time. A correct general form for g^\q,r) in terms of the X 2n , replacing the 
incorrect exp(— ^q 2 X 2 ), is obtained. A simple experimental diagnostic determining when the field 
correlation function gives the mean-square displacement is identified, namely g^(q, r) only reveals 
X 2 if g^\q,T) is a single exponential in r. Contrariwise, if g^'(q, t) is not a single exponential, 
then g^\q,T) depends not only on X 2 but on all higher moments X 2n . Corrections to the crude 
approximation g^(q,r) = exp(— ^q 2 X 2 ) closely resemble the higher spectral cumulants from a 
cumulant expansion of g^(q,r). 



I. INTRODUCTION 



Recently, there has been interest in describing light scattering spectra in terms of mean- 
square displacements of the scattering particles. This description of scattering spectra is 
readily traced back to Berne and Pecora's book 1 , where it is shown that the field correlation 
function g^(q,r) from a quasielastic light scattering spectrum of dilute, true Brownian 
particles is determined by their mean-square displacements ((AR) 2 ) through various times 
r, namely 



where q is the scattering vector. The calculation of Berne and Pecora was meant to describe 
experiments on, for example, dilute solutions of monodisperse polystyrene latex spheres in 
pure water. 

Prior to Berne and Pecora's work, Hallett and students^ used quasielastic light scattering 
spectroscopy (QELSS) to study diffusion of dilute scattering spheres, used as optical probes, 
through complex fluids, namely hyaluronic acid and dextran solutions. The sphere diffusion 
coefficient was combined with the Stokes-Einstein equation 



to infer a viscosity, which was compared with macroscopic measurements. Here ks, R, 
and T are Boltzmann's constant, the probe radius, and the absolute temperature. This 
optical probe diffusion technique was a natural extension of studies by Laurent, et al.— 
of sedimentation of proteins and polystyrene spheres through dextran and other polymer 
solutions. 

The optical probe method has since been extensively been applied to colloid 5 , surfactant^, 
and polymer solution systems^. Techniques for avoiding polymer binding by probes have 
long been understood 9 . Extensive information on the effects of polymer molecular weight 
and concentration on probe diffusion has been obtaine d 10 ^ 1 ; changes in probe diffusion in 
polymer solutions, attendant to changes in solvent quality, have been determined 1 ^. Rigid 
and flexible probes have been compared 1 ^. The diffusion of hard-sphere probes has been 
contrasted with the diffusion of flexible-polymer probes through the same polymer solu- 
tions, solutions being more effective at retarding chain probes than at retarding the motion 
of spheres through the same solutions. 14 Nonexponential spectra have been analyzed by 
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cumulant expansions, Laplace inversion, and multi-mode fitting. The diffusion of polymer 
probes through polymeric matrix solutions, whether observed via QELSS or other tech- 
niques notably fluorescence recovery after photobleaching and pulsed-field-gradient NMR, 
is the topic of polymer self and tracer diffusion, whose literature has recently been system- 
atically reviewed^. 

Since Hallet's invention of the optical probe diffusion method, additional experimental 
methods have become available. QELSS examines single-scattered photons. With increasing 
probe concentration multiple scattering becomes more important, but may be suppressed 
with a two-laser-beam two-detector homodyne coincidence spectroscopy (HCS) apparatus^. 
Popescu, et al.,— have shown how short-coherence- length light can be used to restrict sig- 
nificant scattering to small volumes from which only single scattering contributes to the 
spectrum. Alternatively, diffusing wave spectroscopy (DWS) takes full advantage of the 
properties of light scattered by a profoundly multiply-scattering fluid 18 . All techniques can 
reach the same time scales, because all techniques use the same photodetectors and correla- 
tors, but because DWS responds to the joint motion of many particles it includes relaxations 
due to motions of single particles that are too small to lead to appreciable relaxations in 
QELSS spectra. Note that in the optical probe community 'dynamic light scattering' (DLS) 
refers exclusively to QELSS, while the microrheology community uses 'dynamic light scat- 
tering' as including DWS, so each reference to DLS must be interpreted against experimental 
details. 

Under the cognomen microrheology, an alternative school for the study of optical probe 
diffusion has recently appeared^. Suppose that the mean-square particle displacement of 
particles in an arbitrary liquid were available. This supposition was not advanced by Ref. Q 
in their presentation of eq^ A non-exponential QELSS or other light-scattering spectrum 
S(q, t) could then formally be said to define a time-dependent diffusion coefficient, via 

g?( q ,T)=ex P (-Tq"D(T)), (3) 

where q is the magnitude of the experimental scattering vector and r is the delay time. 
Laplace transform of D(t)t, interpreted as AR 2 , and application of a generalized Stokes- 
Einstein equation in which D and rj are generalized to frequency-dependent forms D(uj) 
and rj(u), might then be proposed as a path to measuring a frequency-dependent complex 
modulus G{oj). 
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This paper treats the field correlation function determined with QELSS. The interpreta- 
tion of experiments using DWS is not considered. Sections II and III of this paper recall at 
minimum length the historical approaches that do lead to eq^under its assumed conditions, 
and that might be said to lead to eq El Section II covers Berne and Pecora's treatment^ of 
the field correlation function for dilute Brownian particles, while Section III covers Doob's 
First Theorem^ for correlation functions of stationary, jointly Gaussian-random, Markoff 
processes. Section IV obtains the actual relationship between the scattering spectrum and 
the mean-square and higher moments of the distribution of particle displacements. If the 
light-scattering spectrum of a suspension of dilute particles is determined by the mean-square 
particle displacement, then it is necessarily also true that the light scattering spectrum is a 
simple exponential. Conversely, if the light scattering spectrum is not a simple exponential, 
for example if it has a secondary slow mode or if it decays as a stretched exponential in 
time at large time, then the spectrum is not determined by the mean-square particle dis- 
placements through different times. Section V analyzes in detail a specific model system, 
thereby clarifying issues related to low- (7 limits and spectral cumulants. A Discussion closes 
the paper. 

II. LIGHT SCATTERING SPECTRUM OF IDEAL BROWNIAN PARTICLES 

This Section considers the light scattering spectrum of a dilute solution of identical true 
Brownian particles, showing that in this special case the QELSS spectrum does determine 
the mean-square particle displacement during different intervals. In a QELSS experiment, 
an equilibrium solution or suspension is illuminated by a narrow pencil of coherent light. 
The light scattered by the solution through a fixed angle is collected. The time-dependent 
intensity I(t) of the collected light is measured. In time domain, the QELSS spectrum S(q, t) 
is obtained as the time correlation function 



Here q is the magnitude of the scattering vector, r is the delay time, and the brackets (• • • ) 
denote a time average over t. Under normal experimental conditions, the time dependence 
of the measured light scattering spectrum is determined by the field correlation function 
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g^'(q, t), namely 

S(q,T) = A[g^(q,T)} 2 + B. (5) 

Here B is the time-independent baseline and A, the amplitude, is a partially apparatus- 
dependent constant. Several normalizations are in common use for A and B. The S(q,t) 
that is directly computed by a real digital correlator is a series of large integers, each of 
which is a sum of large numbers of terms and rij +T being the integer numbers of 

photons counted during short time intervals centered on absolute times j and j + r. 

9 {q-i T ) is i n turn determined by the positions rj(t) of each of the N scattering particles at 
various times. Objects that do not scatter light may affect the spectrum, but only indirectly, 
by modifying the dynamics of the scattering particles. In the limit that the concentration of 
scattering particles is small, the same-time positions of pairs of scatterers are uncorrelated, 
and the field correlation function reduces to 

S (1) (<7^) = (Y,^ex P [ iq .AR l (r)]j. (6) 

Here g^\q,T) is the incoherent structure factor, termed incoherent because it does not 
include terms involving distinct pairs of scatterers. The term af is a normalized scattering 
cross-section for particle i. If all particles scatter equally, af may be absorbed into A. The 
displacement of particle i during r is ARj(r) = r,,(i + r) — Each displacement is an 

integral of the particle velocity v,(s) over times (t, t + r), 

AR,:(r) = | + dsv t (s). (7) 

Further analysis requires a physical description of the particle motions. In the simplest 
idealization, noninteracting probe particles in solution independently perform ideal Brow- 
nian motion as described by the Langevin equation. This special case is treated by Berne 
and Pecorai, who show from the Central Limit Theorem that particle displacements AR, 
over times r much longer than the velocity autocorrelation function's relaxation time, have 
for ideal Brownian particles a probability distribution 



-| -3/2 

G s (AR,r) 



y((AR) 2 ) 



exp [-3(Ai?) 2 /2((AR) 2 )] . (8) 
For Brownian particles, the average displacement over r is 

((AR) 2 ) = 6Dr, (9) 



where D is the single-particle diffusion coefficient of the Stokes-Einstein equation. 
The incoherent structure factor is then 

N 

gM(q,T) = / dAR G s (AR,r)^exp[iq- AR 8 (r)] ~ exp(-Dq 2 r) 
J i=l 



(10) 



An equivalent form is eqQ as seen in Berne and Pecora^. Eq^is sometimes cited as arising 
from Pusey and Tough's treatment 21 of interacting particle systems, but section 4.3.2 of Ref 
L scrupulously footnotes its derivation of eqUJas coming from RefQ 
For simple Brownian particles, the natural logarithm of g^\q,T) is proportional to 
the mean-square displacement of a particle during r, and is linear in r, so the slope 
d log^ 1 ) (q, r)/dr can be used to determine D. In practice, g^ 1 ' is never quite a perfect 
single exponential, and recourse must be had to a more elaborate analysis, such as the 
cumulant expansion 22 . 

As emphasized by Berne and PecoraA, their calculation of g^'(q, t) refers to noninteracting 
particles performing ideal Brownian motion as governed by the simple Langevin equation, 
in which the external random force has a correlation time of zero. This Brownian motion 
idealization is not applicable to any physical system. In particular, the calculation does 
not describe noninteracting spheres diffusing in a small- molecule Newtonian solvent. In a 
real Newtonian liquid, the hydrodynamic wake around a moving particle creates a memory 
effect. The 'random' force on a probe particle thereby gains a long-time memory term, and 
does not have a correlation time of zero. The memory term gives each sphere a velocity 
autocorrelation function that decays at long time as r _3//2 , giving g^(q,r) a tail that does 
not decay exponentially at long times. Paul and Pusey 2 ^ have seen this tail experimentally 
in the Brownian motion of large polystyrene latex spheres. 



III. LIGHT SCATTERING SPECTRUM ON THE BASIS OF DOOB'S FIRST 
THEOREM 

An alternative calculation of the light scattering spectrum follows from Doob's 
treatment^ of Brownian motion and Ornstein-Uhlenbeck processes. Doob's actual interest 
was a particle performing one-dimensional Brownian motion. The particle had a position 
x(t), dependent on time t, and a nominal velocity u(t) = dx(t)/dt. Credible mathematical 
properties of u(t) were shown by Doob to lead automatically to the temporal autocorrelation 
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function (u(t)u(t + r)) . We first recount the key points of Doob's calculation and then show 
how it predicts spectra. 



Restated in modern terms, Doob's First Theorem (Theorem 1.1 of ref. l2Cl) provides 
Let u(t) (with — oo < t < oo allowed) be a family of random variables, with a single 

parameter t labelling the members of the family, determining a stochastic process with the 

following properties: 

(1) The process is temporally homogeneous, i.e., the distribution function for u(t) is 
independent of t. 

(2) If s and s' are two arbitrary, distinct values of t, then u(s) and u(s') have a non- 
singular bivariate Gaussian distribution. 

(3) The process governing the evolution of u(t) is a Markoff process. 

Define m and <Tq by m = (u(t)), o-q = ({u(t) —m) 2 ), and p(t) = ((u(s) —m)(u(s+t)—m)), 
with (. . .) denoting the expectation value. Because the process is stationary, only one value 
for m and one value for o"o exist. Then the stochastic process is of one of two types: 

(A) If s 7^ s', then u(s) and u(s') are independent, Gaussianly-distributed random vari- 
ables, or 

(B) There is a constant (3 > such that if times satisfy si < s 2 < ■ ■ ■ < s n , then 
u(si),u(s2), ■ ■ ■ u(s n ), have an n- variable joint Gaussian distribution, with ((u(si)—m)(u(sj— 
m)) = a 2 exp(-(3 | Si - Sj |). 

To demonstrate this result, Doob observes: We may without changing anything important 
rescale any u(s) so m = and o"o = 1. From (2), the conditional distribution for u(t), given 
a specified u(s), is then 

(27r) 1 /2(i_ p 2)i/2 ex P ^ 2(1 -p 2 ) J { ' 

with p = p(t — s). Consider an ordered series of times t± < t 2 < ■ ■ ■ < £„, and define 
Ui = u(ti), for i e From (2), any pair Ui,uj with i ^ j have the joint distribution 

of eq HU stationarity ensures that p only depends on the time difference | t j — t j \. Eq [TT] 
only gives the two-time distribution function for the Wj. From the Markoff property (3), the 
n-time distribution function is determined, because each is sensitive only to the value of 
the most immediately previous Uj-i, so the Ui have an n-variable joint Gaussian distribution, 
namely 

1 / „,2 m 1 K+ 1 ~ Pi u . 



ex P H/2-ig (i2) 



for pj =p(t j+1 -tj). 

Choosing n — 3, direct integration of eq IT21 shows 



p{h ~ h) = p(h ~ h)p(h - h). 



(13) 



It follows that p is an even function, that p(t) < 1, and therefore either p = (outcome (A)) 
or 



for some positive (3. Eq El is clearly consistent with eq ED Doob implies but does not 
explicitly prove uniqueness. 

It is important to emphasize that it is the Markoff property of u(s) that leads to eqs 
WI\ and [131 The joint Gaussian property of eq [TT] in and of itself does not imply that the 
Ui have an n-variable joint gaussian distribution function; it only implies that any pair 
of the Ui have a two-variable joint distribution function. In a non-Markoff system, u n is 
simultaneously sensitive to all Uj, j < n, not only to the most recent past Uj to be specified. 
The derivation in Section II did invoke the Markoff property of u(s), though not explicitly: 
Equation |H] covertly but necessarily follows from the statement that a displacement over a 
large time can be described as a sum of a series of independent displacements, each over a 
shorter time. The statement that successive Brownian displacements are independent is the 
Markoff assumption, describing the differential pieces of each Brownian displacement; it is 
not the bivariate Gaussian distribution assumption. 

Doob's derivation of his results depends only on mathematical properties (l)-(3), not 
on the physical nature of u(s). The velocity u(s) of a Brownian particle does have these 
three properties, so the Brownian velocity does satisfy (B), but that occurs because the 
Brownian velocity u(s) is stationary, has a two-time joint Gaussian distribution, and is a 
Markoff process, not because u(s) represents the rate at which an ideal Brownian particle is 
changing its position. Correspondingly, any other variable, that has the same mathematical 
properties that u(s) has, will also satisfy Doob's theorem. 

We now consider for which systems the time dependence of the incoherent structure factor 
can be obtained with Doob's First Theorem. The key step is to identify an appropriate 
variable to replace u. An interesting choice of variable is the q th spatial Fourier component 
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of the scatterer concentration 

N 

a <?( s ) = $^ ex p(-q- r i( s )) , ( 15 ) 

the incoherent structure factor being proportional to (a q (s)a*(s + r)). Does a q (s) satisfy 
conditions (l)-(3) of the theorem? 

(1) The a q (s), labelled by the one parameter s, are measured on an equilibrium system, 
so they represent a stationary process. 

(2) At each time s, a q (s) is the sum of a large number N of independent random functions 
exp (— q • Ti(s)). The individual random functions are identically distributed, so by the 
Central Limit Theorem their sum has a Gaussian random distribution. Furthermore, the 
change in a q (s) between two times is determined by 

N 

Aa q (t -s) = a q (t) - a q (s) = ^(exp[-iq • r;(t)])(l - exp[-iq • (r^s) - r^t))]). (16) 

i=i 

The final form is a sum of N terms. Motions of different particles are independent, so 
the N terms are independent from each other. Each term is determined by the particle 
displacement rj(s) — r^t) between times s and t. If all particles are the same, the distribution 
of displacements is the same for every particle, and the same for all particle initial or final 
locations, so Aa q (t — s) is the sum of a large number of independently-distributed random 
variables. From the Central Limit Theorem, Aa q (t — s) therefore has a Gaussian random 
distribution. The convolution of two Gaussian distributions is a Gaussian, so a q (t) = a q (s) + 
Aa q (t — s) also has a Gaussian distribution, and a q (s) and a q (t) have a joint Gaussian 
distribution. 

(3) If an individual particle is performing Brownian motion, and all particles are the same, 
information on the position r^tj) is given by the position at the most recent previous time 
tj-i. If the particles have no memory, once r^tj^) is given, information on the particle's 
position at earlier times rj(i,-_^), £ > 1, gives no further information on how the particle 
moves between tj-i and tj. That is, for Brownian particles the r-j(t) and each of their 
functions is governed by a Markoff process. This Markoff behavior does not follow from 
the Central Limit Theorem rationale used to defend (2), because the displacement of each 
particle between t\ and t 2 could be correlated with the displacement of the same particle 
between t 2 and t 3 . 
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Under these rather restrictive conditions on particle motion, a q (t) has the same mathe- 
matical properties as the u(t) examined by Doob.. In particular a q (ti), a q (t2), and a q (t^) have 
a three- variable joint Gaussian distribution. The autocorrelation function of a q (t) is there- 
fore the same as the autocorrelation function of u(t), leading to an exponentially-decaying 
form 



The incoherent structure factor is almost never a pure exponential. Even the slight 
polydispersity of the spheres in a polystyrene latex preparation leads to a measurable non- 
exponentiality of the incoherent structure factor. Where does this non-exponentiality arise, 
relative to the prediction of Doob's theorem? The key issue is that if the system is poly- 
disperse, then all particles are not the same; some are more mobile than others, and the 
mobile particles continue to be more mobile as time passes. In a series of times ti, t 2 , £3, the 
difference a q (t 2 ) — a q (ti) gives information about the likelihood of changes a q (t 3 ) — a q (t 2 ). 
For example, a particularly small change in a q between times t± and t 2 suggests that the 
fluctuation in a q initially involved an unusually large number of larger, less mobile particles, 
a condition that will tend to persist between times t 2 and £3, so if the change in a q (t) between 
£1 and t 2 is slow, the change in a q between £ 2 and £3 is also more likely than usual to be 
slow. In consequence, for a polydisperse system a q (t) is not governed by a Markoff process. 
Even though its a q (t) at every pair of times has a joint Gaussian distribution, if a q (t) is not 
governed by a Markoff process a q (t) at a trio of times is not described by a three- variable 
joint Gaussian distribution. Without the Markoff condition, an equation like EH but with 
a q (t) replacing u(t) is invalid. Doob's analysis thus does not predict that (a q (s)a*(t)) of a 
polydisperse suspension is a decaying exponential. 

As a further example of a non-Mar koff scattering experiment, consider light scattering 
electrophoresis. In a series of times ti, t 2 , £3, the particle displacement between t\ and t 2 
gives information about the expected displacement between t 2 and £3. In consequence, the 
results here do not lead to the incorrect prediction that the time correlation function in a 
light scattering electrophoresis experiment is a decaying pure exponential. 

The above refers to non-exponential spectra arising from probe polydispersity. However, 
any other effect that leads to a particle motion with memory, such as viscoelastic effects giv- 
ing long-time correlations to the random force on each particle, has the same consequences: 
Memory effects are present, so three a q (t) do not necessarily have a trivariate joint Gaus- 
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sian distribution. The spectrum is not a decaying single exponential. The non-Gaussian 
behavior deduced from this consideration lies in the a q (t), not directly in the distribution 
of particle displacements. However, if successive displacements of every particle had the 
same multivariate joint Gaussian distribution function, then three a q (t) would also have a 
trivariate joint Gaussian distribution, eqEHwould follow, and the spectrum would be a pure 
decaying exponential. Therefore, if the incoherent structure factor is nonexponential, not 
only trios of a q (t) but also trios of successive particle positions do not have three- variable 
joint- Gaussian distribution. 



IV. LIGHT SCATTERING SPECTRUM BY DIRECT CALCULATION 

Superficially, it would appear that the light scattering spectrum could be calculated via 
a Taylor series expansion of eqEl namely 

t=l \n=0 

n being the Taylor expansion index, with the N particles potentially not all being the 
same. Inversion symmetry causes all terms odd in ARj(r) to average to zero. For identical 
Brownian particles that follow eq |H1 and m a positive integer 

<(q • ARi(r)) 2m ) = q 2m 2 m {Dt) m {2m - l)(2m - 3) . . . (1). (19) 

On substituting eqE3into eqEH the factors (2m — l)(2m — 3) . . . cancel the odd factors of 
n!, while the term 2 m is cancelled by the factors of 2 in the even factors of n\, leading to 

JV oo , _ 2-t\n N 

9 {1) M = J2J2- z ^t l ^J2^-^/ 2 ^ (2°) 

8=1 n=0 1=1 

where X 2 = ((Axi(r)) 2 ) and where Axj(r) = q ■ ARj(r). For dilute Brownian particles, the 
Taylor expansion approach recovers the results of the previous sections. 

What if the diffusing particles do not execute simple Brownian motion? The difference 
between the actual g^(q, r) and its Mar koff-pro cess component is seen by the factorization 

g ( 1 )(q ) r ) = exp(-g 2 X2/2) (^exp (iq ■ AR(r)) + q 2 X^ '/2^j } . (21) 
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Within the average, expanding separately the exponentials in AR(r) and X 2 , sorting terms 
by their order in q, and re-exponentiating, one finds the approximation 



9 {1) (q,r) 



exp 



-tx* 



{1 + ^(((Ax^)) 4 ) - 3(X 2 ) 2 ) - 0(q 6 ) + ...} (22) 



Define X n = ((Ax(r)) n ) . The terms 1 + aq 4 are reorganized by multiplying by unity 
in the form exp(ag 4 ) exp(— aq 4 ), factoring an exp(+ag 4 ) out of the term in braces, and 
simultaneously multiplying within the braces by the Taylor series for exp(— aq A ), following 
which the g 4 term vanishes from within the braces. The lead correction to the exponentials 
is then q 6 ( — 2X 2 ^ + X 2 X 4 — 2X 6 )/48. Iterating the process of factoring out lead terms as 
exponentials, eq 12*01 can be rewritten through 0(q 8 ) as 

/ V2 Q~V2 2 \ finv? 3 



cxp 



2 X 2 4 (X 4 



24 



3X 2 ) + ^ (30X 2 - 15X 2 X 4 + X 6 ) 



720 



-q 



,(630X 5 



420X 22 X 4 + 35X 42 



28X 2 X 6 - X 8 ) 



, 1 + 0(q w ) + ... (23) 
40320 1 

Equation [2S1 is the general expansion for the incoherent structure factor in terms of 

moments of the distribution of particle displacements, complete through order g 8 . The 

mean-square displacement by itself determines only the lead term of g^(q, r). 



For nearly Brownian particles, X 2n increases approximately as t n . If g^(q,r) is not a 
simple exponential, the terms of eqESlof order q 4 and higher do not vanish, and \og(g^ (q, r)) 
does not give the mean-square particle displacement. It might superficially appear that com- 
plications arising from terms of 0(q A ) and higher could be avoided by moving to sufficiently 
small q. Ignoring the modest detail that in conventional units under any credible experi- 
mental circumstances g > 1, the superficial appearances are misleading. To clarify this, we 
examine a special case. 



V. SPECTRUM OF A BIDISPERSE SYSTEM 

Suppose that the field correlation function is exactly 

g {1 \q,r) = A ie x P (-D iq 2 T) + A 2 ex P (-D 2 q 2 r). (24) 

Here A± and A 2 are two mode amplitudes and D\ and D 2 are two diffusion coefficients. 
This g^(q, r) would arise from a bidisperse suspension of Brownian particles; it would also 
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arise from a monodisperse suspension in which Brownian particles found themselves in two 
different environments. Each mode separately corresponds to particles whose motions are 
describes by a Gaussian Markoff process with X 2 = 2DiT. 

Equation |^] is a uniform function of q 2 r. If one changes the scattering angle, thereby 
changing q, and simultaneously changes the units of time, so that q 2 r does not change 
(equivalently, if one plots the spectrum against q 2 r rather than against r), then the shape of 
the spectrum is independent of q. By going to lower scattering angle, one reduces q, causing 
each of the two decay modes to appear at later times, but other than a change by a constant 
multiplicative factor in the temporal positions of all spectral features, the spectral lineshape 
does not change when q is changed. 

The field correlation function is numerically well-behaved so it and its logarithm have 
Taylor series in q 2 . The simple Taylor series is 

The exponential of the Taylor series of the logarithm of eqEUis 

(i)/ \ /a A s ( (MDi + A 2 D 2 ) 2 A 1 A 2 {D 1 -D 2 )\ 2 2 \ f . 
g w {q,T) = (A 1 +A 2 )esxp[ q t+ (q r) +... , (26) 



v (Ai + A 2 ) * 2(A 1 + A 2 ) 2 7 

Note that the g 4 terms of eqs 1213 and [2H1 are not the same. Even in our special case, in which 
the motion of each particle is separately totally characterized by its mean-square displace- 
ment as a function of time, the scattering spectrum includes nontrivial terms determined 
by the difference (AxA 2 {D\ — D 2 ) 2 )/(2(Ai + A 2 ) 2 )(q 2 r) 2 between the actual mean-fourth 
displacement, and the mean- fourth displacement expected for Gaussian particles having the 
same mean- square displacement. 

The series in eq recovers the original function and all of its properties. In particular, 
g^(q,r) as plotted against q 2 r is invariant under a change in q and, therefore, egH3] as 
plotted against q 2 r is invariant under a change in q. The relative importance of the q 2 and 
high-order terms does not change as q is changed, so eq EHl does not tend toward a single 
exponential in r as q — > 0. If q is changed by a factor /, the time scale on which eq HE\ 
relaxes changes by a factor f~ 2 , but the spectral lineshape does not change. 

In our special case the A{ and the Di are simply numbers, not functions of time. Therefore, 
comparison with the orthodox spectral cumulant expansion 

9 <%r)=exp(f: £tp2), (27) 
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where the Kj are the spectral cumulants, is directly possible. The series in eqs ESI and 1771 
(after taking logarithms) are both convergent power series in r, so they are equal term by 
term, leading to 

K = \og(A 1 + A 2 ) (28) 

The first cumulant, which determines the initial decay of the field correlation function, is 
a weighted average of both diffusion coefficients; the initial decay rate does not give the 
diffusion coefficient for the faster mode. The second cumulant is the q A term. Except 
in systems containing highly monodisperse particles in very simple fluids, higher spectral 
cumulants are not negligible. Because the spectral lineshape is invariant to a change in q 
the higher cumulants cannot be made less important by going to smaller q. 

Equations EH refer to the model given by eq UJ] In general, the cumulants may be 
obtained from logarithmic derivatives of g^(q, r), namely for n > 1 

-Qj) log(s (1) (g,r)). (31) 
In terms of the particle displacements, the first few cumulants are 

K 1 = lim (£!f!_£lMI^) + ...| (32 ) 

1 i-o I 2 dt 24 dt 1 

and 

Kl = lm (l S ^-l^- r i + (33) 
t-o I 2 dt 2 24 dt 2 1 

For the Brownian particles of the special case, X 2 ~ t 1 and X 4 ~ t 2 , so only the X 2 term 
contributes to K\ and only the X 4 term contributes to K 2 . This simplification cannot 
be made without some specification of the physical properties of the system. Also, in 
interpreting the limit t — > 0, the original field correlation function g^(q, r) is only measured 
at a series of times r that are much larger than the relaxation time of the diffusing 
particles, so the limit also only reaches to times ^> r#. Care must then be taken to interpret 
the derivatives in the above equations as examining correlations between particle positions at 
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some time t and particle velocities at later times t + r; without this interpretation anomalous 
outcomes may follow 2 ^ 

How does taking the q — > limit affect eq 1261.'' Suppose the limit of small q is taken. In 
this case, a point with fixed r = r refers to a point earlier and earlier in the spectrum's 
decay. If q is truly small, g^(q, r) = 1 to within the precision of the experiment. At slightly 
larger q, the decay of g^(q,r) will from to r G be linear in r, and will therefore within 
the precision of the experiment be indistinguishable from a pure exponential having K\, 
above, as its relaxation rate. Note, however, that this description of g^'(q, t) as a purely 
g 2 -dependent exponential only refers to very early times. At later times, the 0(q ) and 
higher terms of g^'(q, t), whose coefficients are independent of r and are therefore present 
at all times, rise above the limits of experimental precision. No matter how small the non- 
zero q, if t is increased enough that the field correlation function relaxes to zero, all the 
higher-order-in q terms of eqEHl become substantial. 

VI. DISCUSSION 

The objective of this short paper was to show how the incoherent structure factor g^ (q, r) 
measured by light scattering spectroscopy is related to the distribution of particle displace- 
ments during the interval r. To put the discussion into its correct historical context, the 
standard result, eq Q for ideal Brownian particles, was obtained both from the Langevin 
equation and from Doob's Theorem. For ideal Brownian particles — \og(g ( - 1 \q, r)) is deter- 
mined by the mean-square particle displacements ((AR) 2 ). 

As seen in eqESl lead terms for the general relationship between g^'(q, r) and the average 
particle displacements were obtained. Except at very small times, g^(q,r) is in general 
determined not only by the mean-square particle displacements, but also by the higher 
moments of the particle displacement distribution. Plots of — log(g^(q, r)) against t, taken 
to long times, do not reveal the mean-square particle displacement through long time, except 
in the special case of ideal Brownian particles, for which g^(q,r) is a pure exponential. 
However, if the spectrum is an exponential, there is no need to determine a limiting slope 
at large t, because the slope is the same at all times, and therefore the short-time behavior 
determines X 2 correctly. There are non-exponential QELSS spectra of probes, in complex 
fluids, that have already been analyzed by applying eq^ It is important to emphasize that 
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nothing is wrong with the underlying spectra. Re-analysis of those spectra using orthodox 
QELSS approache d 22 ' 24 may well yield interesting information about the fluids. 

There remains the case of particles whose dynamics are complex at short time, but whose 
motion when adequately coarse-grained appears approximately Brownian. This case would 
arise, for example, for Brownian particles passing through an ordered or random potential 
energy field having a well-defined longest length scale £. Particle motions over distances 
^> £ would smooth over the variations in the potential energy field, so that if trapping was 
not an issue the particle motions over very long distances would be approximately Brownian 
and eqQ would again be applicable, even though the short-time motions were not Brownian. 

However, there is no guarantee that the dynamics of a complex fluid are appropriately 
characterized as having a longest length scale, or that the spectrum is a single exponential 
at long times. Many complex fluids have hydrodynamic interactions, whose longest range is 
given by the Oseen tensor with its r _1 interaction, leading to relaxation effects on all distance 
scales. Correspondingly, relaxations are found with all time scales, leading to a g^\q, r) 
having a long-time stretched-exponential-in-time relaxation, exactly as seen experimentally^ 4 - 
in some cases. It is, of course, possible to plot the log of a stretched-exponential against 
time out to long times, but that plot shows a smooth curve that lacks a long-time straight 
line limit. One could always fit the last few points of such a to a straight line, but that 
line is only the local tangent of a smooth curve that continues to bend at longer times. The 
slope of that fitted tangent line is determined by the largest r at which one measured g^\ 
and is therefore purely an artifact of the fitting process. 

In the above, it has been shown that in general QELSS spectra do not determine mean- 
square particle displacements X 2 . The special-case exception to this rule, in which QELSS 
spectra do determine X 2 , refers to identical particles performing Brownian motion in a 
simple homogenous fluid. In this special case, the spectrum is a pure exponential, and all 
spectral cumulants other than Kq and K\ vanish. 

How does this result affect interpretation of the literature? The optical probe diffusion 

nnnnn 

literature, e.g., refs. I8U9I )13U14I)25L appears to be unperturbed by the above results, because 
in the optical probe literature eq. ^is not used to interpret QELSS spectra. Instead, QELSS 
spectra are fit to exponentials, to cumulant series or sums of cumulant series, or to special 
functions, or are subject to Laplace inversion. With any of these fits, the objective is to 
generate a short list of parameters that characterize the spectra, and to examine how the 



16 



fitting parameters are related to other solution properties or to fundamental models for 
polymer solution dynamics. Because eq^was not used, it is not significant that it is not 
correct for probes in complex fluids. 

A rather different situation is found in the microrheology literature, that is experiments 
tracing themselves back to Ref. lij in which eq ^ is systematically applied to interpret 
light scattering by optical probes. Much of the microrheology literature uses diffusing wave 
spectroscopy; the above analysis did not determine if eq^is correct for DWS measurements. 
However, some of the microrheology literature, e.g., Mason, et al.— proposes that QELSS 
'can also be used to measure the mean square displacement of probe particles', and therefore 
can be used in the same way that DWS can be used, a proposal that as seen above is incorrect. 
It is important to emphasize that the difficulty is in applying eq Q to the spectra, that 
nothing is wrong with the spectra themselves, and therefore that spectra can be correctly 
reinterpreted using methods developed for studies of optical probe diffusion. As examples 
of the approach envisioned by Mason, et al.— , consider: 

1) Dasgupta, et alM, studied probe diffusion in polyethylene oxide: water using DWS, 
QELSS, and eq[T] Their QELSS data, using very large (0.97 /iin) probes in 4wt% polyethy- 
lene oxide- water, was reduced via eq^to a time- dependent mean-square displacement, whose 
angular dependence scales as q 2 . As noted by Dasgupta, et al., at short times their QELSS 
and DWS data begin to differ, and their QELSS data was not reported for times shorter than 
a large (by QELSS standards) 10 mS. The deviation is what would plausibly be observed if 
DWS followed eq[T]even though QELSS follows eqEHl Our own 29 QELSS data on probes 
in aqueous 900 kDa polyethylene oxide, extended down to a few /iS, confirms that optical 
probe spectra in these solutions at short times are radically nonexponential, and therefore 
that eq^is not applicable for probes in polyethylene oxide:water as studied by QELSS. 

2) van der Gucht, et al.— used QELSS to study diffusion of 250 nm probes in solutions of 
the associating polymer bis(ethylhexylurido)toluene as a function of polymer concentration. 
Their probe spectra resemble those of probes in solutions of hydroxypropylcellulose^, probes 
showing a near-exponential relaxation at small polymer concentration and an additional, 
prominent, long-lived mode at elevated polymer concentrations. Here g^(q, r) is interpreted 
in terms of X 2 , with early and late regions of the spectrum interpreted in terms of diffusive 
decays (i.e., X 2 ~ t 1 ). 

3) Kang, et al.— report on the diffusion of spherical probes through solutions of fd- viruses, 
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which are effectively rigid rods, using fluorescence correlation spectroscopy, video microscopy, 
and QELSS to measure diffusion of optical probes in different time regimes. They assume 
that QELSS gives the mean-square particle displacement, an assumption based on an ex- 
pansion in terms of small wave vectors, and assert that they are able to enter this regime 
experimentally. In this regime, as discussed above, g^ l \q-,T) is necessarily indistinguishable 
from a single exponential that has decayed little from its average value. 

Finally, in searching the literature on this topic, note that the optical probe and microrhe- 
ology literatures are appreciably non-communicating; a citation search of either literature 
will largely not couple into the other. 
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